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Abstract 

We study the effects of noise and decoherence for a double-potential 
well system, suitable for the fabrication of qubits and quantum logic el- 
ements. A random noise term is added to the hamiltonian, the resulting 
wavefunction found numerically and the density matrix obtained by av- 
eraging over noise signals. Analytic solutions using the two-state model 
are obtained and found to be generally in agreement with the numerical 
calculations. 

In particular, a simple formula for the decoherence rate in terms of the 
noise parameters in the two-state model is reviewed and verified for the 
full simulation with the multi-level system. The formalism is extended 
to describe multiple sources of noise or different "dephasing" axes at the 
same time. Furthermore, the old formula for the "Turing- Watched Pot" 
effect is generalized to the case where the environmental interactions do 
not conserve the "quality" in question. 

Various forms for the noise signal are investigated. An interesting 
result is the importance of the noise power at low frequency. If it vanishes 
there is, in leading order, no decoherence. This is verified in a numerical 
simulation where two apparently similar noise signals, but differing in the 
power at zero frequency, give strikingly different decoherence effects. A 
short discussion of situations dominated by low frequency noise is given. 

1 Introduction 

In the study of macroscopic quantum mechanics and quantum logic devices in 
particular, the question of decoherence and its effects remains an important, if 
not the most important, issue. Inevitable disturbances from the environment of 
the quantum device, and correlations established with variables of the environ- 



merit, will limit the length of time for which the system under study acts as an 
isolated quantum system. 

In previous work we studied, by means of numerical simulation, how deco- 
herence affects and interacts with the behavior of some quantum logic devices 
PP. The decoherence was modeled by a simple random noise signal presumed 
to act on a certain component of the system - in the SQUID corresponding to 
flux noise. In the present work we wish to consider more generally the connec- 
tion between decoherence, the type of random noise signal, and the parts of the 
system upon which it acts. We reach a number of interesting conclusions on the 
effects of these factors. 

A further interest of such studies is its connection with the "quantum mea- 
surement problem" , and we are able to illustrate quantitatively phenomena like 
the "Collapse of the Wavefunction" and the "Watched Pot Effect" through the 
simulations. 

2 Hamiltonian 

We shall examine the hamiltonian 



which is suitable for representing a "qubit" and through manipulations on x ext , 
the quantum logic operations NOT and CNOT pQ. 

The hrst part of the expression represents a double potential well problem 
illustrated by Fig [T] The Hj\f(t) represents a small, time dependent, random 
noise term. This will be used to simulate the decoherence effects. These could 
be due to either actual noise sources in the laboratory like stray fields, or true 
quantum mechanical decoherence due to unavoidable interactions of the simple 
idealized system with other variables. 

The parameter x ext controls the asymmetry of the two potential wells, x exf = 
being the completely symmetric situation where the energy splitting between 
the two lowest levels ui tunnel is only due to the tunneling through the barrier. 
The height and width of the barrier are controlled by ft. In the previous work 
it was found that the values ^ = V Q = 16.3; ft = 1.19 give behavior suitable for 
quantum logic elements and we shall use these values in our simulations here. As 
seen on the figure, with these values there is a well defined pair of levels below 
the barrier. At the same time the barrier is low, giving a relatively large tunnel 
splitting. This permits the adiabatic sweeps, on which the gate operations are 
based, to be fast. For these values one has w ttmT!e z = 0.0044. 

For the SQUID, these parameters are related to the inductance L and capac- 



itance C by u = V w 1030 y L /^ H ■ The parameter ft characterizes the critical 

current I c and is given by ft = 2v £ c . (We use K, c = 1 units.) The overall 

energy scale is given by E a = 1 / V LC and the time unit by 1 /E a ~ \J LC. Thus 
the time unit is 1.0 x 10~ 12 'seconds x (L/pH C/pF) 1 / 2 . With L = 400pH and 



H = —-— + V Q {\[{x-x ext f]+ftcosx}+H N {t) , 



(1) 
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Figure 1: Double well potentials (V term) of the hamiltonian Eq[TJ with the first 
four energy levels indicated. The left panel shows the symmetric configuration 
with x ext = while the right one shows an asymmetric configuration with 
x ext _ 0.0020. The behavior of the two lowest levels is well represented by an 
effective hamiltonian H = ^cV, where V x is given by the tunneling energy and 
V z by the asymmetry of the wells. 

C = Q.lpF for example, Lo tun nei — 0.0044 corresponds to ~ 4 x 10~ 7 eF. The 
time unit is then 6.3 x 10~ 12 sec. [2]. 

The numerical methods are as outlined in [TJ, as introduced in |3]. A large 
basis of harmonic oscillator states is employed, and the resulting sparse hamil- 
tonian matrix is inverted to find the energy eigenstates. The time evolution is 
found by repeated application of (1 + ^iHAt)/(l — ^iHAt). 

3 Two- State System 

The previous studies have shown that there is a range of parameter values where 
the two lowest levels for EqTTJcan be adequately represented as an effective two 
level system with hamiltonian H = ^crV. 

We use an (x,y,z) coordinate system for the various vectors in the following 
way. V is approximately in the 'z' or 'up' direction when x ext is relatively large 
so that the potential is asymmetric and the energy level splitting is essentially 
determined by the distance between the bottoms of the potential wells. V is in 
the 'x' direction or 'horizontal' when x ext — and the energy splitting is deter- 
mined only by the tunneling through the barrier. Specifically, one determines 
P] the components of V to be V z « 2V x ext x c and V x = uJtunnei, the tunneling 
energy. The quantity x c is the value of the x-coordinate where the wavefunction 
is centered when the wavefunction is localized in just one of the potential wells. 
For the parameter values used here this is always close to 1, about 0.9. The value 
of UJtunnei, namely =0.0044, is obtained from the numerical solutions as the level 
splitting for x ext = 0, with the chosen parameters fi = V Q = 16.3; j3 = 1.19. 

According to [T] , deviations from the two-state model begin to appear when 
the level splitting for the lowest pair becomes comparable to the distance to the 
next set of levels. Various tests showed that the effective two-state description 
was good up to about x ext « 0.01, where the level splitting was 0.30. 
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To represent the noise in this approximate two-state picture, we assume 
an additional term ^crB in the two-state hamiltonian with a time dependent 
random field B. The effective total two-state hamiltonian is then 



H = ±<tV+ ±<tB. 



(2) 



B thus represents a rapid small stochastic addition [5] to the V. Just as for 
the components of V, it is possible to establish relations between the B and the 
parameters of Hj^. These will be given in section [7] 

4 Density Matrix 

Our object of study is the density matrix of the system. In our random noise 
model we generate it by inserting a particular realization Af a (t) of a noise signal 
in the hamiltonian EqQ] and then evolving an initial wavefunction to a final 
time with this hamiltonian H a . Repeating the process, and averaging over 
many different noise realizations, N, gives a final density matrix: 



Although the evolution is unitary for a given J\f a (t), the averaging process is 
not and the procedure evidently violates unitarity. The resulting density matrix 
will thus exhibit a time-dependent decoherence, which we wish to examine. 

In particular we would like to study the decoherence for the effective two- 
level system comprised of the lowest states of the hamiltonian, with its 2 x 2 
density matrix. However, for the calculations with the full EqQ] it must be noted 
that the resulting numerical density matrix Eq[3]is a function of two continuous 
variables in p(x',x). Thus in principle Eq|3] refers not only to the two lowest 
states of the hamiltonian, but to the entire Hilbert space of the continuous x 
variable. 

Nevertheless, in the previous studies it was found that in the region of pa- 
rameter space where the two-state description is valid, the Hilbert space of 
lowest states is accurately spanned by one and the same set of wavefunctions 
even as the parameters of the hamiltonian are varied and the energy eigenstates 
change [4 . Therefore wavefunctions beginning in the two-state space remain 
there, and the results of evolving with EqQ] can be expressed in the two-state 
space by projecting them onto the corresponding pair of wavefunctions. This 
is the procedure we will follow in the numerical evaluation of the 2x2 density 
matrix. 

On the other hand, for the simplified two-state or spin 1/2 system repre- 
sented by Eq[2l the density matrix can be predicted analytically, as will be 
explained in section [5j This analytical result may be compared with that for 
the numerical procedure just described, where we insert the noise signal in EqQ] 
To carry out this comparison one needs a translation between the noise signal 
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inserted in H^f and the equivalent B in Eq[5] This translation will be given in 
section [7] 

The two-state density matrix can be parameterized in terms of Pauli matrices 

as 

p=§(l + P-<7). (4) 

With our coordinate conventions P is in the "up" or "down" or ±z direction 
when the wavefunction is localized in one of the potential wells of FigfTJ and it 
is in the "horizontal" or ± x direction when the wavefunction is the symmetric 
or antisymmetric combination of states in both potential wells. |P| = P = 1 
corresponds to pure state and P < 1 a mixed state. 

In the two-state model the "polarization" vector P evolves according to [14] 

P = PxV L>P T , (5) 

where the — D Pt term represents the quantum damping or loss of coherence. 
D is the decoherence rate and D~ Y is the decoherence time. This term reduces 
the components P transverse to the direction chosen by the noise or decoherence 
interaction. 

Although in the present paper we concentrate on situations with constant 
parameters in the hamiltonian (aside from the noise term), in previous work time 
dependent V 's were used. These may be employed to realize adiabatic quantum 
logic gates [T] , and via the adiabatic method to give a direct measurement of D 
by a turning on-and-off of the classical-quantum transition [8] . 

5 D in the two-state model 

We now consider how the decoherence or damping and the associated parameter 
D are generated by the random noise in the two-state picture of Eq[2j The V 
are assumed constant or slowly varying over the time period necessary to find 
the decoherence effects induced by the random fields. The method is essentially 
the same as in Ref p] but we would like to present it in a somewhat more general 
form. This is nicely done by following Ref [9]. One considers the density matrix 
p=i(l + P-cr)in the Heisenberg representation. The P are also assumed to be 
slowly varying so that at a given time one may assume a certain initial density 
matrix which we call p(0). 

In Eq[3]one may denote each of the terms on the right as p a , so that p = ~p a . 
We then examine the density matrix equation p = —i[H(t),p] term by term 

p° = -i[H°(t),pP], (6) 

iteratively: p{t) « p(0) + p\ + pi- The different orders refer to how many factors 
of the small B occur and one compares equations with the same number of -B's 
[5] . Since we will assume the average value of the random fields B to be set to 
zero (or incorporated in the main hamiltonian), the first order contribution to 
the average p is zero: pi = — i[erB, p(0)] = 0. We thus consider the second order 
for p 
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p a 2 (t) = -i±[<7B ,tf] = -\[*B a , [aC a ,p(0)]] (7) 

where we introduce the shorthand C(i) = J Q dt'B(t') and use pi = — [<rC, p(0)] . 

Now in performing the average over the random fields in Eq[71 we will have 
to do with bilinear correlators between the vector components of B, C. These 
can be thought of as tensors of the type XjYj , where the average is over different 
realizations of the random fields. With just one random field in a fixed direction 
e, one has 

X®Y=Ce®e t (8) 

for the tensor, where C represents some quadratic correlator. 

By making repeated use of [ctX,ctY] = i2aX x Y and comparing the 
resulting coefficients of cr, Eq[7] implies 

P = BC(ex(exP)). (9) 

To evaluate the correlator BC we use the basic result from the theory of sta- 
tionary random noise that the quantity (J Q B(t')dt') 2 = 2tA, where A is the 
time integral of the autocorrelation function 



A = I B(0)B(t)dt . (fO) 

Then 



BC = B(t) J B{t')dt' = |^ B(t')dt'^j =A. (II) 

As is evident when taking the scalar product with e, the cross product 
expression represents the component of P perpendicular to the respective e, i.e. 
with the component along e removed. So we may call ((eP)e — P)= — Pt with 
"T" for transverse. 

We therefore obtain that the D in Eq[S]is given by 

D = A (12) 

E q[T2l has the typical form of a dissipative parameter related to the integral over 
the autocorrelation function for a fluctuating quantity [10 . 

One can also handle more than one noise source by this method. With two 
independent random fields, along the fixed directions ei and e2, we will have a 
tensor of the form 

X (g) Y = dex <g> ei f + C 2 e 2 ® e 2 f . (13) 

That is, different noise effects can be added up separately and do not interfere 
with each other. This is due to the assumed independence of the two noise 
signals. Note the e x and e 2 need not be orthogonal. 



1 This is the first point where one uses the properties of the a as Pauli matrices. Otherwise 
they could have been any group generators, in which case the cross product formulas are 
replaced by bilinear products of structure constants [9] 
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Thus according to Eq[T51 for the case of two independent noise signals Eq[5] 
generalizes to 

P = [Ui((eiP)ex - P) + -A 2 ((e 2 P)e 2 - P))] . (14) 

Calling (eiP)ei P = — Pxi,and (e2P)e2 P = Pt2 we thus find that 
the old [13] rotation- decoherence equation Eq[5]now becomes 

P = PxV-AP T1 -D 2 P T2 , (15) 

with 

D x = Ai D 2 = A 2 (16) 

Concerning dimensions, B is an energy in view of its role in the hamiltonian, 
and so A via EqUUlis also an energy, which matches the units of D. 

The conclusion Eq |15l could also have been arrived by an intuitive argument 
where one imagines turning off one noise signal for a short time and then repeat- 
ing the process with the other noise signal turned off. One would then obtain 
first the results of Eq[S] with respect to one direction and then with respect to 
the other direction. The net result would correspond to E q[T5l 

As said, the e's can be in any direction, not necessarily orthogonal; and since 
the independent addition of signals in Eq[T3l holds for any number of different 
signals, evidently E q|15l can be extended to more than two noise signals by 
simply adding a — D Pt term for each. It is of course essential that the different 
noise signals be uncorrelated. 

The length of P is always decreasing with an equation of the type Eq|T51 and 
these considerations can be given an interpretation in terms of the increase of 
entropy [15]. In particular we note a formula for the decrease of the length of 
P (squared): 

1 1 P 2 = P P = P ( Dl P T1 + D 2 P T2 ) = -APti 2 ~ A>P T2 2 , (17) 
z at 

for the case of two independent decoherence or noise signals as in E q[T5l 

Concerning the necessity of considerating more than one noise source, it is 
possible or indeed likely that there are noise and decoherence effects on more 
than one aspect or component of a system. For the SQUID, for example, there 
might be external flux noise giving fluctuations in x ext and at the same time 
effects in the Josephson circuit amounting to fluctuations in j3. Such a situation 
would be described by E q[T5l 

With only one noise "axis" the decoherence problem may be viewed as one 
of random phases or rotations around that axis, as was done in [T| for flux noise 
in the SQUID. This is sometimes called "dephasing", and D is the associated 
diffusion parameter. Thus it may be said that Eq[T5l gives the formula for the 
evolution when there is more than one "dephasing" axis. 
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6 Noise Models 



We shall consider different types of noise signals 7V(t). Noise signals are charac- 
terized by their power spectrum V(lo), which is the fourier transform of the au- 
tocorrelation function J\f(0)J\f(t), namely V{oS) = J* °° dtcoswtN(Q)N(t). One 
such power spectrum is "white noise" with a cutoff ui c : 

This is a constant or white noise spectrum up to frequency w c , and then falls off 
as 1 /w 2 at high frequency. This was used in a simplified way in pQ . In section lTTI 
below we will discuss some examples using this kind of noise. 

In addition, we also use planckian spectra corresponding to 1, 2, or 3 spatial 
dimensions. With these we can to try to represent a thermal background and 
as will be explained below, allows us to conveniently compare cases with and 
without noise power at zero frequency. They also allow a simple variation of the 
frequency content by changing T which is a "temperature " parameter charac- 
terizing the frequency spectrum, with an exponential cutoff above u) ~ T. For 
n=l we take the planckian form 

V{.) = ^- i n=l. (19) 

We have fixed the normalizations in E q[T9l and E q[T8l such that V(0) = 1. 
This condition corresponds to setting the time integral of the autocorrelation 
function equal to one: 

V{0)= dtAf(0)jV(t) = — S a / dt J\f a (0)J\f a (t) = 1 , (20) 

Jo ^ J o 

where N a (t) is a particular realization "a" of the signal, and we have generated 
N such signals. This convention is particularly convenient since according to 
Eq ll2l the decoherence parameter D is simply proportional to just this time 
integral of the autocorrelation function. With this standardized normalization 
of the noise signal we will regulate its overall strength in different applications 
by an adjustable small coupling factor. 

The planckian forms corresponding to spatial dimensions n=2 and n=3 can- 
not be normalized in this way since they contain extra powers of uj, giving 
vanishing noise power at lj = 0. We therefore choose to normalize them such 
that the variance of the signals, which is a measure of their magnitude, is equal 
to that for n=l at the same T. 

This leads to the following definitions 

- 5 ^\ n = 3. (21) 
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The variance of the signal corresponds to the integral over the power spec- 
trum, so with these normalizations the integral J Q duiV(ui) is the same for all 
three signals in EafT^l and EoJ5TJ 

Specifically, one finds 



For T=0.025, which we shall use in many of the simulations, this gives a 
typical signal strength of yJM 2 = 0.16. As we shall see below in section [TU1 the 
case of noise signals with the property f dtAf(0)J\f(t) — 0, as in Eql2~Tl are of 
special interest. 

To numerically generate a particular signal J\f a (t) as a realization of one 
of the above models, we work with a discretized fourier space where Af(t) = 
£ w A/Le lw *. A given random set a of the AC then determines the signal. We 
have used two procedure to generate the J\f u . 

In the random phase procedure, one takes ^JV(ui) times a random phase e 
with 9 chosen randomly and uniformly in the interval > 8 > 2n. 

In the gaussian modulus method one chooses a real and imaginary part of 
Af u randomly from a gaussian distribution with zero mean and variance V(uj). 
Both methods lead to an ensemble with the power spectrum ~ V(u>). We find 
that both procedures lead to very similar results in our simulations. 

Figure shows some samples of the different noise signals. The upper curve 
is for Eq[P8l One sees the presence of the high frequency tail in the power 
spectrum. The middle curve is for the planckian n=3, Eql2"T| and the lower curve 
for the planckian n=l, Eq[T9l These latter two differ in that for n=3 there is no 
power at zero frequency, while for n=l it is present. The temperature T=0.025 
or uj c is the same for all three. It will be seen that the amplitude of the signals 
is roughly in accord with Eq[22l with y/J\f 2 = 0.16. 

7 Relation Between Hj\f and B 

In this section we determine the relations between the random B of the two-state 
description Eq[2] and the noise term in the full hamiltonian Eq[T] Given these 
relations, we will be able to compare the results of a full numerical simulation 
with Eq[T] and the analytical predictions from EqrT2"l 

As explained below, noise in the 'z' direction arises through fluctuations in 
x ext , and in the 'x' direction, fluctuations in j3. In the case of the SQUID these 
correspond to noise/decoherence in the external flux and in the junction circuit, 
respectively. It is quite plausible that they should correspond to independent 
noise sources. 

7.1 Noise in the "z-direction" 

We thus consider adding a noise term representing a fluctuation of x ext in Eq[TJ 




(22) 



H M = 5x ext V a x 



8riM{t)V x. 



(23) 
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Figure 2: Samples of noise signals Af(t). Upper curve: Cut-off white noise from 
Eq llSl with cutoff oj c = 0.025. Middle curve: planckian noise for n=3 as in Eql2~T1 
with T=0.025. Lower curve: planckian noise for n=l as in Eq[T9l with T=0.025. 
The signals are normalized to have the same variance or integrated noise power 
but differ in having zero (middle) or non-zero (upper and lower) noise power at 
zero frequency. The high frequency tail present in the power spectrum of the 
upper, non-planckian, curve is evident. With L = AOOpH and C — O.lpF for 
the SQUID, the time unit in this and all further plots would be 6.3 x 10~ 12 sec. 

We call the fluctuation Sx ext and for a particular noise signal " a" it is represented 
by SijM a (t) with 8rj an overall coupling strength. By using the identification 
x ~ x c a z (see section [3]) connecting the operator x of the full hamiltonian with 
the a of the spin 1/2 picture, on sees that this Ujv is equivalent [5] to the B 
field of the two-state model as : 

B z = 2V x c 5r]M , (24) 

with other B components zero. We therefore predict, according to Eq[T2l that 
when using a noise signal normalized with 'P(O) = 1 as explained in section [6l a 
decoherence constant from the full simulation with EqQ] 

POO 

D= B z (0)B z (t)dt = 4:(V x c Sj]) 2 = 0.89 x 10 3 {8-qf , (25) 
Jo 

where in the last step we insert Vq — 16.3, x c = 0.9. 
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7.2 Noise in the "x-direction" 



Switching to the "x-direction" , we consider the effect of small variations of (3 in 
EqflJ amounting to fluctuations in the barrier potential. Since /? regulates the 
tunneling, a change in /3 induces a change in the tunneling energy. Thus in the 
two-state hamiltonian Eqf5]a variation in ft induces (apart from an unimportant 
constant shift) a change in the a x term, a B x . Other components of B are zero. 

To obtain the magnitude of B x [T2], one notes that the coefficient of <j x is 
the tunnel splitting; V x + B x = w tU nnei + 8cj t unnel- Thus B x = Sui tunne i = 



z 




Figure 3: Typical configuration of vectors used in the simulations. With no 
damping the polarization vector P rotates around V , so with V in the x- 
direction there are oscillating components P z ,P y . For D not too large these are 
damped oscillations as the component of P perpendicular to the noise direction 
JV or B is reduced. For very strong D, P is "pinned" along the noise vector, 
corresponding to the "Collapse of the Wavefunction" . The initial P is of unit 
length and along the z-axis for Figs[5] and [6l while for Fig[7] it lies at 45 degrees 
in the x-z plane. 

Although there are analytic methods of evaluating tunneling energies, they 
involve exponentially sensitive effects and it seems best to find d " t -" ei from 
our previous numerical results. In the vicinity of /? = 1.19, (and with Vb = (a = 
16.3 where we carry out our simulations) examination of Table II of [1] yields 
^gasi ~ _o.l8 near /3 = 1.19. Thus, for fluctuations around j3 = 1.19 the B 
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field associated with a change in /3 is 



B x = 8(3 ™' = (-0.18)5/3 = {-Q.\%)Sp o N , (26) 

We write the fluctuations in (3 as S(3 = 8(3 N ', with N one of our noise signals 
and J/3q a strength parameter. According to Eq[T2l when Af has the integral of 
the autocorrelation function normalized to one, 

D = 0.032 5 ft . (27) 

For the numerical calculations, we wish to insert the same 8(3 fluctuations 
in Ecjl] Since x 2 /2 is the leading term in cosx which contributes to the level 
splitting, a small variation in (3 can be adequately represented by adding a term 
8(3V x 2 /2 in the full hamiltonian. We thus should set, in EqQ] 

H N = -8/3 V \x 2 = -8p a N V \x 2 . (28) 



7.3 Comparison of x and z-noise 

Our two examples of 'z-noise' and 'x-noise' are instructive in that the effective 
B fields for the two-state system E q|2~il and Eq|2l)l and the associated D's, Eq|2*5l 
and Eq |27l differ widely for similar noise signals. The difference originates in 
the fact that the 'z-noise' operator is linear ~ x, while the 'x-noise' operator is 
quadratic ~ x 2 . The two lowest states comprising our two-state system have, 
with x ext = 0, odd or even parity. Thus they will be strongly connected by the 
parity odd x operator. For the x 2 operator on the other hand the change in the 
energy splitting is only due to the different shifts of the two states under x 2 and 
so is much smaller. For x ext ^ but small the situation will be qualitatively 
similar. 

Thus the symmetry properties associated with the noise or decoherence can 
be important. With specific devices this may suggest a choice of operating 
conditions where the effects can be minimized. 



8 Tests by Numerical Simulation 

We turn to numerical tests of the formalism and the above results Eql25l and 
Eql2"7lof the spin-^ picture. We check them against results from the full hamil- 
tonian Eq[T]with the corresponding Hj\f, namely Eq|2Ulor E q|2"51 inserted in the 
hamiltonian. We find the decoherence parameter D from a fit to the damping of 
|P| in the numerical calculations and compare it with the analytic predictions. 
For our various tests we chose a coupling 8rj such that the decoherence time 
1/D is on the order of some thousands of time units. For the typical SQUID 
parameters of mentioned in the introduction this corresponds to a decoherence 
time in the tens of nanoseconds. 
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Figure 4: Excitation to higher states. When the noise power spectrum contains 
frequencies above that of the principal level spacing, here 0.4. the diagonal 
elements of the 2x2 density matrix fall below 1/2. When it does not contain 
these frequencies the relaxation is to 1/2. 



8.1 Parameter Ranges 

There will be certain conditions on the range of parameters where the analytical 
arguments can be applied. 

One condition concerns large perturbations, where the nature of the system 
is changed substantially. These are evidently beyond the scope of the theoretical 
method, and would have to be dealt with simply by simulation alone. This is 
particularly relevant in our example of /3 fluctuations or 'x- noise' where the 
tunneling energy is very small and j3 fluctuations have a relatively large effect. 
This is discussed in sect 18.41 

A second condition concerns the frequency spectrum of the noise, as char- 
acterized by the T or the w c of the noise power spectra. If the dominant noise 
frequencies are too low, on the order of, or less than the important frequencies 
of the system, this is in contradiction to our assumption that the system can 
be taken as constant over many cycles of the noise [T3] . In sect[TT] we shall 
discuss some features of low frequency noise. On the other hand, if the frequen- 
cies in the noise are too high, one risks exciting higher states and so violating 
our assumption of an effective two-state system. Since the important system 
frequency we deal with is uitunnei — 0.0044, and the distance to the next set 
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V 


Pinitial 


D 


(0.0044,0.0,0.0) 


(1,0,0) 


0.00134 


(0.0044,0.0,0.0) 


(0,0,1) 


0.00122 


(0.0044,0.0,0.006) 


(1,0.0,-1)72 


0.00118 


(0.0044,0.0,0.02) 


(1,0,0) 


0.00141 



Table 1: Study of the constancy of D for various starting conditions. A n=l 
planckian noise in the 'z-direction' with T = 0.05 is used. The coupling is 
5r) = 0.00128, for which the prediction Eq[25]is D = 0.00146. In the first two 
entries V is along the x-axis while the starting P is first along the x-axis and 
in the second line along the z-axis. In the third entry V is close to 45 degrees 
with the z-axis and P starts orthogonal to it in the lower hemisphere. In the 
last and fourth entry V makes an angle of only 12 degrees with the z-axis and 
P starts along the x-axis. One observes that at about the 20% level the value 
of D is the same for the various examples and near to the predicted value. The 
notation for the vectors refers to {x, y, z). 

of levels is about 0.4, one has the condition 0.0044 << T,uj c « 0.4. This 
criterion applies of course for our typical configurations and would have to be 
adjusted for others. 

The excitation to higher states above that of our effective two- level system 
can be detected in that the diagonal elements of the 2x2 density matrix will relax 
not to 1/2, but to some smaller value, reflecting the missing probability. We can 
perform an amusing numerical experiment on this by using a 'rectangular' noise 
power spectrum V(u)) which is exactly zero above some to = T. The resulting 
behavior of the density matrix is shown in Fig |4] for various values of T. 

If the noise power does not contain frequencies equal or higher than those 
needed to go to the next set of states, the density matrix should relax to 1/2. 
With our typical parameters, this frequency is about 0.4. In Fig|4]it is seen that 
when the maximum noise frequency is below this value the relaxation is indeed 
to 1/2. With higher frequencies present the relaxation is to less than 1/2. (That 
it falls to 1/4 reflects the fact that the numerical calculations were carried out 
in a 4-state basis.) 

Although our various noise power spectra do not have a simple rectangular 
cutoff, we anticipate a similar behavior with our exponentially cutoff planckian 
spectra. We do indeed find with low T's and couplings that are not too strong 
that calculations performed with a two-state basis give the same results as basis 
with more states, indicating that the excitation of higher states is not important. 
The results we present below are always under these conditions. 

8.2 State Independence of D 

Probably the most striking aspect of the theory is that, with the given as- 
sumptions, the effect of the environment, here represented by the noise, can be 
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V 


T 


D 


Pn 


(0.0044,0.0,0.0) 


0.0125 


0.0011 


0.50 




0.025 


0.00128 


0.50 




0.05 


0.00132 


0.50 




0.10 


0.00145 


0.50 




0.20 


0.00157 


0.45 



Table 2: Temperature dependence of the decoherence parameter D determined 
from the numerical simulations using noise signals in the 'z- direction' with 
the power spectrum E q[T9l and coupling 8-q = 0.00128. According to E q[25l D 
should be temperature independent and have the value 0.00146. The final value 
of the density matrix element pu should be ^ if there is no excitation of higher 
states above the effective two-state system. The last entry for p n shows some 
deviation from this as T approaches the principal level spacing 0.4. 

represented by a single parameter D. Most importantly, D does not depend on 
the state of the system. This means that Eq[5] is a linear equation for P . This 
would not be the case, for example, if one had D = D(P). Linearity for the 
evolution of the density matrix is a general result of quantum mechanics and it 
is non-trivial that our approximate method respects this. 

We can provide some tests of this feature in the full numerical simulation. 
One may take a given noise signal with a given coupling and determine D for 
different configurations of the potentials and of the initial state. In TablcQ] 
we show some examples, using 'z-noise' with the parameter Srj — 0.00128, for 
which, according to E q|2"5l the prediction is D=0. 00146. 

The first two entries are for symmetric potentials, x ext = 0, that is with 

V in the x-direction. For the first line the system is started with P also in the 
x-direction, that is with a wavefunction which is a linear combination of equal 
wavepackets in the left and right potential wells. In the second line the system 
is started with P in the z-direction, corresponding to the initial wavefunction 
concentrated in only one potential well. In the next two entries we show the 
case of some asymmetric potentials, so that V lies at some angle in the x-z 
plane. In the third line V is at a 45 degree angle, and in the fourth line close 
to the z-axis at a 12 degree angle. 

8.3 Test for 'z-noise' 

We now examine the nature of the evolution with "noise along the z-axis" . With 

V along the x-axis, Fig[3] shows this situation with P in some general direction. 
In a first test we start P along the z-direction again with coupling 5rj = 

0.00128. One anticipates that P will rotate around the x-axis while its length 
decreases. Inserting Eql23lin Eq[[J evaluating the density matrix and extracting 
P one finds the results shown in Fig[5] The left panel shows the projection of 
P on the z-axis, with the expected damped oscillations. The right panel shows 
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the decrease of the total length of P 
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Figure 5: Simulation with Eq[T] using "z- noise" due to x ext fluctuations as given 
by Eql23l The configuration corresponds to Fig[3]with P started along the z- 
axis. Left: Projection of P on the z-axis showing damped oscillations. Right: 
Damping of the total length of P . A fit for the length with e~ Dt yields D = 
0.00132, while the prediction from E q[25l is D — 0.00146. The error bars indicate 
the dispersion of values expected for an average made from N samples, namely 

\J Y>i (P a — P) 2 , where P a is the quantity being plotted, resulting from a given 
noise signal 'a'. 

The value of D emerging from the simulation is found by fitting e~ Dt to the 
curve of the right panel, yielding D=0. 00132. On the other hand the prediction 
from Eql25l is D=0. 00146. Thus there is good agreement. 

These plots were produced using the n=l planckian noise Eq[T9l It is a 
feature of the prediction Eq |25l with fixed 5rj that the D resulting from the 
simulation should be given by V(0) and independent of T, even though a change 
of T implies a change in the frequencies present in the noise. 

To test this, Table [2] shows the values of D resulting from various T's used 
in simulations with n=l planckian noise. The first and last values of the tem- 
perature are near the margins of the allowed region 0.0044 << T,u c << 0.4. 
One observes approximate agreement with the analytic prediction D=0. 00146. 
At the highest temperature the excitation to states above the lowest two states 
manifests itself in the relaxation of the density matrix element to a value less 
than J- 

These results are meant simply as a check on our numerical and mathemat- 
ical methods and should not be taken to mean that the decoherence is temper- 
ature independent. To the contrary, a temperature dependence is rather to be 
expected. With the normalizations we have adopted, it would be represented 
by a temperature dependence of the coupling parameter 5rj. 



8.4 Test for 

Fluctuations in j3 correspond to changes in the tunneling energy, a a x operator 
in our notation. As was discussed in section [7731 fluctuations in j3 have relatively 
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small effects compared to similar fluctuations in the x ext parameter, leading to 
a much smaller D for (3 fluctuations or 'x-noise'. 

If for the purpose of comparing the analytic and numerical predictions, we 
were to increase the coupling constant 8/3 in the simulations, there is the dif- 
ficulty that using a larger 8(3 will induce substantial changes in the small tun- 
neling energy, violating our assumption that the fluctuations do not essentially 
change the system. For example, we work with V x = uj tun nei = 0.0044, and 
according to E q|2rJl the B x induced by a 8(3 is (— 0.18)<5/3 o N. If we wish to have 
B x « 0.0044 with the typical size of noise signals 0.16 (see sectO, then we 
require 8(3 « 0.01, which implies D on the 10 -6 level, according to EaBTl 

While such small effects can be unstable and difficult to detect numerically, 
we show in Fig[6] a run for x ext = with the coupling 8(3 a — 0.01. The tem- 
perature for n=l planckian noise was set to T=0.05. A fit gives D — 5 x 10~ 6 , 
while the prediction of Eq[27] is D = 3 x 10~ 6 . Although there is qualitative 
agreement, it is evidently not so precise as in the case of the z-noise. This may 
be due to our various approximations as well as numerical uncertainties. 




Time Time 



Figure 6: Relaxation due to (3 fluctuations or "x-noise". Left: Projection of 
P on the z-axis. Right: Damping of the total length of P . For the configuration 
of Fig[3]with P started along the z- axis. A fit to the length with e~ Dt yields 
D = 5 x 10~ 6 , while the prediction from Eq|25l is£> = 3xl0~ 6 . A weak 
coupling S(3 = 0.01 implying weak damping must be used to avoid a large shift 

8oJtunnel /^tunnel ■ 



9 "Collapse" and "Watched Pot" 

Since the beginnings of quantum mechanics the "Collapse of the Wavefunction" 
has been assumed to be an instantaneous and somewhat mysterious process, not 
describable in a physical manner. However, as was noted in the beginning of 
studies of decoherence [2], the action of the external influences, here simulated 
by the noise, may be quantitatively examined and the induction of a "Collapse 
of the Wavefunction" -like behavior of the system studied. This occurs through 
the effect of the —DPt term, which tends to eliminate off- diagonal matrix 
elements of the density matrix and so to reduce the density matrix to a clas- 
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sical probability distribution, an incoherent mixture of the states of a definite 
"quality" , the eigenstates of a z . 

As the external influence, represented by D, becomes very strong, the "Col- 
lapse of the Wavefunction" occurs very rapidly. In addition, the further evo- 
lution of the system is inhibited and in the limit can be essentially stopped 
completely 1 14) -much like the 'fixing' of the state by a 'measurement'. This is 
the "Turing" or "Watched Pot" effect. In terms of Eq[5]this occurs because the 
large D "pins" P to the noise axis so there is a very small Pt and very little 
subsequent damping. 

9.1 Strong Damping with jV, V Perpendicular 

With our present system we can illustrate and study such effects quantitatively. 
Let the initial P be chosen to have length one (pure state) and be oriented at 
an angle such that both P z and P x components are substantial. Fig|3] indicates 
such a starting condition. If there is a strong noise along the z-axis, P x will 
rapidly disappear: there is a "Collapse of the Wavefunction" to an incoherent 
mixture of "up" and "down" states. After this "collapse" the large value of D z 
inhibits further evolution- the "Watched Pot" effect. 

FigOshows a simulation for such conditions. The length of P is plotted. One 
observes a rapid decrease of the length until only the original P z survives. The 
wavefunction "collapses" in a very short time . After this, the length remains 
"frozen" to this value-the "Watched Pot" effect. Increasing D gives a faster 
"collapse" and a stronger "watching". As with the other simulations, the full 
hamiltonian Eq[T]is used here. 

The parameter characterizing "strong" or "weak" damping is the ratio D/w, 
where uj — V is the energy splitting of the two levels or the magnitude of V . 
This follows from Eq[5l where dividing through by V yields a dimensionless time 
parameter tu> and the parameters D/w. For small D/w one has weakly damped 
precession of the P vector around V. For D/u> >> 1 one enters the "Watched 
Pot" regime. 

Examination of the solutions of Eq[5] for V in the x- direction with V x = 
^tunnel and P initially in the z- direction leads to the behavior |14) 

P z « e-W) (29) 

for P initially along the z-axis. This result follows from combining the z and y 
components of Eq[5]to give the equation P z + DP Z + V X 2 P Z =0 and solving for 
large D [TT] . For very large D the graph in Fig[7] will decrease extremely slowly 
(see Fig 6 of pQ). 

9.2 Strong Damping with J\f, V Nearly Parallel 

In the above case Af or B and V were at right angles. This is the usual "Watched 
Pot" configuration as it has been studied for chiral molecules, neutrinos, or 
the SQUID with a symmetric potential and dominant flux noise [14] pQ. In 
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Figure 7: The "Collapse of the Wavefunction" followed by the "Watched Pot" 
effect. The plot shows the evolution of the length of P for the configuration 
of Fig[3] with the initial P of length one and at 45 degrees in the x-z plane. 
The P x component "collapses" to zero quickly, giving an incoherent mixture 
of "up" and "down" states. This is then followed by a near "freezing" of the 
evolution. These effects result from the strong noise (or "measurement"), here 
with D/oj « 100. Note the short time scale, the time unit is the same as in the 
other plots, namely with L = AQOpH and C = O.lpF for a SQUID, the unit is 
6.3 x 10" 12 sec. 



all these cases one has- or it is assumed-, that the observable "quality" in 
question (chirality, flavor, flux state) is conserved in the interaction with the 
environment. Since we take this 'quality' to be represented by an eigenstate of 
tr z , B is necessarily in the z-direction. 

However, one may also envision the opposite situation where Af and V are 
in the same direction. This would be the case, for example, with a strongly 
asymmetric potential for the SQUID [x ext ^ 0) and dominant flux noise. (This 
appears to be the case for the experiment mentioned in ref |13j.) 

In this case the analog to Eql29l for the long term behavior of the parallel 
component after the perpendicular components have been "collapsed" is found 
to be 2 

Pparallel « P^raUel^^ ■ (30) 

■^parallel ^ s the component remaining after the initial "collapse" and V 2 the 
square of V . This result follows from writing Eq[5] as a set of linear equa- 
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tions, taking P ~ e * and solving the determinental equation for the smallest 
eigenvalue A. 

For D >> V one recovers E q|2"51 For D << V there is a new regime with 
a suppression due to the large value of V, with the behavior ~ e~ tDe where 
& = V x /V. Naturally in the limit where V and Af are exactly parallel, P will 
remain constant, as one sees setting V x = in EqIBTil 

10 Vanishing Noise Power at Low Frequency 

In the numerical work so far we have assumed that the time integral of the 
autocorrelation function, or equivalently 'P(O), is non-zero. However a very 
interesting special case occurs when it is in fact zero as for Eq[2~Tl Then Eq[T6l 
implies 

/>oo 

D = / B(0)B(t) = (31) 

J o 

To check if Eq l31l holds with the full hamiltonian we perform simulations 
comparing noise signals with planckian power spectra for n=l and n=3. If 
E q|3"T1 indeed holds there should be little or no relaxation in the n=3 case. Since 
our different signals have been normalized in E q|19l and E q|21l to have the same 
variance N 2 , and the square root of the variance is the typical size of the signal, 
our comparison is with signals of the same 'strength'. 

Employing these signals for noise in the z direction in Eq l23l we obtain the 
results shown in Fig[S] The configuration is as in Fig[3] with P started in the 
z-direction and the parameters the same as for Fig[5j The left panel of Fig[8] is 
for n=l noise and shows the expected relaxation of P z . The right panel is for 
n=3 and shows no apparent decoherence or damping effect in the oscillations of 
P.- 

An examination of the length of P shows that while it goes to zero for n=l, 
for n=3 it declines by only about 1% over the length of the run. 

Thus despite the at first sight similar appearance of the lower two signals 
in Fig[2] there is a great difference in the decoherence they induce and the ex- 
pectation from Eql3T1 is indeed verified in full calculations with Eq[TJ The exact 
vanishing of D in E q|31l is presumably a consequence of the use of second order 
perturbation theory in the derivation and should be lifted in higher order. A 
similar run for n=2 shows in fact some weak decoherence, with D at the ~ 10~ 5 
level, perhaps indicative of the higher order effects. 

11 Very Low Frequency Noise 

Eq[5] gives a local-in-time description in terms of a single "frictional" or "dissi- 
pative" constant D. D does not depend on the state of the system (the value 
of P) or on the internal hamiltonian (the value of V). D depends only on 
the properties of the noise, and our simulations have quantitatively verified 
this. This simple structure has its origin in the "random walk" behavior of 
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Figure 8: Identical runs for P z but with two kinds of noise. Left: Planckian 
noise power spectrum with n=l (lower curve of Fig [2|, having non-zero power 
at zero frequency. Right: Planckian noise power spectrum with n=3, having 
zero power at zero frequency (middle curve of Fig[2]) . The two noise signals are 
adjusted so as to have the same variance or integrated noise power. There is a 
striking difference in the damping or decoherence, even though the equality of 
variances assures signals of the same "strength" . The conditions are the same 
as for Fig[5j 



the integral of the noise signal J* B(t')dt' . If B is a random variable then its 
integral represents a random walk and, as was used in arriving at E q|10[ one 

has (J * B{t')dt') 2 = 2tA = 2tD, with A the time integral of the autocorrelation 
function of the noise. This is the result of a random walk with steps given by 
V2A. Therefore f*B(t')dt' is gaussian distributed with a variance increasing 
in time as ~ t, which finally leads to the exponential in time P = e~ Dt . 

If, however, the noise is dominated by very low frequencies [13] such a local- 
in-time description is not relevant and the methods presented in the above 
discussions should not be expected to apply. For times t small compared to 
the time scale on which the noise B varies, J Q B(t')dt' is not undergoing a 
random walk. Indeed, in the limit a particular realization B a (t) may be taken as 
effectively constant and one has f Q B a (t')dt' — tB a (0) instead of the i 1 / 2 random 
walk. The decoherence will now result from the average over the different B a (0). 

A simplest case to analyze in this way is that of subsection 17.21 where the 
noise vector B is parallel to V (junction circuit noise for the SQUID). V gives 
the rotation rate of P in the z-y plane for a P started along the z-axis. Thus 
F->y + fi with random B gives a spread of the orientation of the P increasing 
in time and so a decreasing length of the averaged P or decoherence according 
to 



P = cosJ B(f)df « cosB(0)t = J cosBt Prob(B)dB t«t nolse (32) 

where t no i se is the time scale over which B changes and Prob(B) is the prob- 
ability distribution for the B. With a gaussian Prob(B) one obtains P = 
e -( 1 / 2 )B 2 t j n g enera i f or sma ll times one will expect P « 1 — \B 2 t 2 . 
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One sees that when t « t no i se the relevant quantity is the magnitude of 
the fluctuations of B, namely B 2 . As is evident, for example, in the discussion 
around E q|2"2l this is not the same as the V(u) = 0) entering into D. Rather B 2 
is given by the integral over 'P(w), Eql22"1 

In Fig|9] we show the result of some simulations with very low frequency 
x-noise compard to a high frequency case. 




5000 10000 15000 20000 

Time 

Figure 9: Simulations with low and high frequency x-noise compared. The 
upper curve shows the effect of a "rectangular" noise power spectrum with 
cutoff ui c = 2 x 10~ 4 and the lower with uj c = 1 x 10~ 2 . The low frequency case 
exhibits the expected P ~ e — const. t behavior anc j the high frequency case the 
P ~ e -const.t behavior of our main discussion. The dashed and dotted lines are 
fits for the two types of behavior. 

We use a cut-off white noise spectrum similar to Eq[T51 However to make 
the situation with the noise frequencies completely clear we use a "rectangular" 
spectrum with V = 1 up to u) c , where it is then abruptly cut off. The typical 
time t no i se characterizing a change in a noise signal is then less than or on the 
order of u) c . The upper curve is for the small w c = 2 x 10~ 4 or t no i se > 5000. In 
order to make the effects evident the relatively large value Sn = 0.1 was used. 
The behavior P ~ e -const.t - 1S geen anc j a (dashed line) gives 2.0 x 10~ 8 for 
the constant. Using Eql26lto find the effective B, and the integral over the power 
spectrum to find the variance, one predicts const. = (1/2)B 2 = 2.9 x 10 -8 , in 
reasonable agreement. The lower curve, on the other hand, is with the higher 
frequency noise ui c = 1 x 10~ 2 = 1/100. One observes the simple exponential 
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behavior (dotted line), with the coefficient D — 4.1 x 10~ 4 in approximate 
agreement with formula Eq[25l for D which gives D = 3.2 x 10~ 4 . That the low 
frequency case has the smaller effective decoherence in this example is due to 
the fact that with the "rectangular" noise spectrum and V fixed to one below 
u> c , the variance B 2 decreases as oj c decreases (Eq l22"|) . While for the purposes 
of comparison we have kept V at w = constant for the two cases, the behavior 
under experimental conditions will of course depend on the actual amplitude of 
the low frequency noise. 

We now turn to the perhaps more experimentally relevant situation of sub- 
section im , as represented by Fig[3] While the V is still along the x-axis, the 
noise is now along the z-axis, ("flux noise" in the SQUID ). P is again started 
along the z-axis. 

The noise perturbation along the z-axis leads to two differences with the 
x-noise case just discussed. One is that the shifts in the rotation frequency are 
now quadratic in B. One has V — > \JV 2 + B 2 s» V+ Thus the magnitude 
of the alteration of V is less than in the x-noise case, but will always have the 
same sign, leading to an increase of the average frequency of rotation, by -^yB 2 . 
This average change is not relevant to the decoherence, but its fluctuations, 
characterized by B 4 — (B 2 ) 2 , will play the role of the B fluctuations characterized 
by B 2 in the x-noise case. 

The second difference has to do with the fact that since B is no longer 
parallel to V , the plane in which P rotates is tilted somewhat with respect 
to the y-z plane. This will lead to some oscillatory effects in the projections of 
P , but not to an increasing-with-time decoherence arising from the spread in 
rotation frequencies. 

We thus again concentrate our attention on the spread in rotation frequen- 
cies. We repeat the arguments around Eq[32l with the role of B(t) now played 
by B 2 (t)/2V. EqGH now becomes 



P = cos( ~^r)dt' « cos(-t) = J cos(-t) Prob(B)dB (33) 
= cos(±arctan(at)) ^ /4 , 

where we have taken a gaussian Prob(B) and a — B 2 /V. For small t one has 
Pal - (3/8)(B 2 /V) 2 t 2 . As before, B 2 is the relevant parameter of the noise, 
but now it enters as (B 2 ) 2 since we are now concerned with the fluctuations of 
B 2 instead of those of B. 

In FigfTUl we show z- noise simulations for P with a low, w c = 8x 10~ 5 and a 
higher u c — 8 x 10~ 4 frequency rectangular noise spectrum. The Sr] parameter 
of Eql24"lis 2.5 x 10 -3 . One notes the above-mentioned oscillations. The envelope 
to the upper curve (dashed curve) can be fit with Eql33l using a — 5.3 x 10~ D , 
while the prediction with the parameters used is a — 6.4 x 1CP 5 . Checks on the 
dependence of a with respect to the magnitude and frequency of the noise signal 
are in agreement with the theory. The lower curve is fit by a simple exponential. 
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Figure 10: Simulations with a very low oj c = 8 x 10 -5 , and higher lu c = 8 x 10~ 4 , 
frequency z-noise. The dashed line fitting the upper curve uses Eq l33l with 
a = 5.3 x 10 -5 . It exhibits the 1 — const.t 2 behavior at small times, with 
const. ~ (B 2 ) 2 . The lower curve is fit with a simple exponential. 



The fact that the decoherence here depends on B 2 or the integral over V(uj) 
leads to amusing consequences in the n=2,3 "planckian noise" cases, where the 
V(0) determining D is zero. As the temperature is lowered and B 2 becomes the 
relevant quantity, the decoherence could increase. We find this is indeed the case 
in simulations, with the decoherence increasing substantially for temperature 
few x 10~ 3 , to above the D ~ 10~ 4 level. 

A further interesting point to be noted is that in the general formalism the 
length of P is always entropically decreasing (see Eq [T7| : while here, as in the 
oscillations in FigfTOl P can in fact increase. 

In general concerning very low frequency noise, we can conclude that there 
will be a spread of oscillation frequencies or "dephasing" at times t « t no i se , 
leading to a decoherence which is related to the variance of the noise signal B 2 . 
as opposed to the integral over the autocorrelation function determining D. It 
appears in different ways according to the type of noise, as in our examples of 
x- or z- noise. While with high frequency noise the t << t no i se regime may 
be unobservably small, with very low frequency noise it could dominate the 
relaxation of P, as for the upper curve of Fig|9l 

In this section we have retained the assumption of a small noise signal: 
B/V << 1 and treated a constant V. In the low frequency limit for the noise 
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it would also be possible to treat a large B, and a slowly moving V, using the 
adiabatic approximation, but this is beyond the scope of the present discussion. 

12 Conclusions 

In summary, we have reviewed and extended the formulas relating random noise 
signals and the decoherence parameter for the effective two level system formed 
by the lowest levels of a double potential well. The random noise is applied 
as a fluctuation of various parts of the full multi-level hamiltonian. Numerical 
simulations verify the theoretical predictions, in which the relevant characteristic 
of the noise is the power at zero frequency. In particular different kinds of noise 
signals are applied and we find reasonable agreement between the simulations 
and the analytic predictions. 

This agreement includes the size of the decoherence parameter, as well as its 
independence on the state of the system and on the frequency composition of 
the noise. Effects related to the excitation of levels beyond those constituting 
the two qubit-like states are discussed. The results confirm the possibility of 
describing the dynamics with an equation with a single decoherence constant D. 
"Collapse of the Wavefunction" behavior as well as that of the "Watched Pot 
Effect" are induced for large values of D. In general one finds that, within the 
allowable parameter ranges, the full complexity of the SU2 structure involving 
the noise vectors, the V and the P arises from the two-state reduction of the 
full hamiltonian. A particularly striking result, in agreement with the theory, 
is the small size or absence of decoherence when the noise signal has no power 
at zero frequency. 

The noise power at zero frequency appears in fundamental treatments [TB] 
[TT] of the dissipation problem as the level of white noise. As is familiar from 
Johnson noise, for example, the existence of a classical dissipative parameter 
like resistance is associated with a constant noise spectrum at low frequency. 
We thus arrive at the perhaps obvious interpretation that to have a small D 
such classical dissipative effects or parameters should be kept to a minimum. 

An interesting general question concerning our approach is to what extent 
classical random noise can be used as an adequate simulation of decoherence. 
A real, classical signal, that is an unquantized [18 field, causes as many transi- 
tions 'up' as 'down' between any pair of states. This means that "spontaneous 
emission" or "relaxation" , where, even in the absence of any noise, an upper 
level emits a quantum (e.g. a phonon, or in vacuum a photon) and transits 
to the lower state is absent. On the other hand, the neglect of spontaneous 
emission may not be important in practice. Quantum logic devices should have 
short time scales of operation, while spontaneous emission rates will be slow. 
Study of this question would require a detailed analysis for specific devices. 

While our main results are in the context of standard high frequency noise, 
we also offer a short discussion and simulations of noise signals dominated by 
low frequencies. It is found and verified by the simulations that in this case the 
relevant parameter becomes the variance or integral over the power spectrum 
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of the noise. 
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